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Abstract 

Within the framework of the Fock-Feynman-Schwinger (equivalently, worldline) 
casting of QCD a four-point process, mediated by a closed quark loop, is studied in the 
Sudakov kinematical region, first perturbatively and subsequently with the inclusion 
of non-perturbative contributions through the use of the stochastic vacuum model. 
Deformations of the, first order, perturbative results are determined. Finally, the 
resummation of the leading logarithmic terms is performed and physical implications 
of the result are discussed. 
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1. Introductory remarks 

Sudakov behavior was originally formulated in connection with a vertex function in QED 
[1]. It pertains to a situation where a very large momentum transfer is imparted on the 
fermionic line by a photon and accounts for the suppression of probability that photon 
radiative emission is not detected above a given (infrared) scale. One of its significant impli- 
cations for pQCD pertains to exclusive processes and, in particular, to the extent by which 
perturbative factorization holds for that corner of phase space where the momentum of the 
hadron is basically carried by a single constituent quark (so called Feynman picture). This 
corresponds, in effect, to a low-x situation, given that it implies the vanishing of the mo- 
mentum fraction carried by the other, accompanying, valence quark(s). In the perturbative 
context the aformentioned infrared scale should, of course, be possibly close, but certainly 
above Aqcd- In practice, it corresponds to the (inverse) separation scale between the par- 
ticipating and the low-x quark(s). The central objective of this work is to study possible 
modifications induced by non-perturbative contributions to the result of Refs. 2-4. 

The methodological tools for our analysis are provided by: a) The Fock [5] -Feynman [6]- 
Schwingwer [7] (FFS), also known as worldline, casting of QCD [8-10], on account of certain 
advantages it offers in connection with the application of eikonal approximation methods in 
QCD [11] and b) the stochastic vacuum model (SVM) [12-16], which has been constructed for 
the purpose of confronting nonperturbative issues in QCD. This scheme is founded on three 
basic axioms whose aim is to define the vacuum structure of the theory, several consequences 
of which have been succesfully checked against lattice results [17]. 

In this work we address the Sudakov suppression issue in an 'dealized' context which 
pertains to the four-point process depicted in Fig. 1. Once isolating the kinematical corner 
which ensures Sudakov kinematics for the process, namely large momentum transfer at each 
one of the four vertices, we proceed to study the problem concerning the extent to which 
non-perturbative contributions affect the relevant, Sudakov result. 

Following a brief exposition of the FFS- worldline formalism regarding relevant basic tools, 
we shall proceed, in Section 3, to determine the Sudakov kinematical region, as it pertains to 
the four-point process under discussion and study its perturbative implications. Section 4 is 
devoted to the quantitative study of lowest order, non-perturbative contributions to the four- 
point process (always in the Sudakov kinematical region) and determine the modifications 
induced by them. Once this is done a resummation of the overall soft contributions will be 
carried out (section 5), which will facilitate direct comparisons with the pQCD result [2]. 
Concluding remarks will be presented in the last section. Finally, the technical manipulations 
involved in the proof of a set of equations in the main text are relegated to an Appendix. 

2. General formalism 

Consider the n-point Green's function entering a QCD process defined by a product of 
local currents J^(x) = ip (x) T^ip (x) , where T M is an appropriate element of the Dirac-Clifford 
algebra. One writes 

^>n-/W i X n) ' ' ~, X l) = {J^i{ x n) ' ' ' Jy. n {%l))if),A 

— Tr c (iG(x 1 ,x n \A)r iXn iG(x n ,x n ^ 1 \A)r iXn _ 1 ■ ■ ■ iG(x 2 ,x 1 \A)T lxl ) A . (1) 
In the FFS-worldline formulation the propagators entering the above relation are given by 
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[8-10] 



iG(x,x'\A) = J™dTe- m2T J Vx(t)e- l *f° dti2{t) [m-^-x(T)} 



x<f>{x] [1/2] Pexp 



x(0)=x',x(T)=x 



ig 



dtx ■ A \ . 



o 



(2) 



In the above relation cp[V2] enters as a factor which accounts for the spin of the propagating 
particle modes. More specifically, spin effects are represented by the term exp (j J T dt o^F^ , 
which should be positioned inside the path ordered Wilson exponential. It is possible, on 
the other hand, to factorize the spin effects via a partial integration [17] and embody them 
on the so-called spin factor [18] which is given by 



$[x]^ = Pexp r-£dta^„[xU , 



(3) 



where = ^[x^(t)x u (t)—x u (t)x^(t)} is a 'torsion' tensor and are the spin-1/2 generators 
for the Lorentz group 1 . In the context of the eikonal approximation where spin effects are 
absent, it can be set to unity. On the other hand its presence will become important in 
connection with the non-perturbative considerations that will eventually enter our analysis, 
at which point it will be explicitly taken into account. 

Under suitable parametrizations the four-point Green's function, to which we specialize 
our considerations from hereon, can be cast in the following form (note that the Sj, i = 1, 4 
serve to parametrize the, closed, quark contour -see Fig. 1- with T = s 4 being the 'total 
time') 



G 



fj,4 ■■■m 



■i=4 J /„■>_ ,rr,_ 



X 



S i*-imi[z] (Tr c Pexp 



z(0)=z(T)=x 
T 

ig ' 



Jo 



dtz ■ A 



where 



i=4 A 



(4) 



(5) 



Going to momentun representation one obtains, after taking into consideration momen- 
tum conservation around the fermionic loop, the following expression for the Green's function 



G> 4 ... W (P4, • • -,Pi) 



i=l 



z(0)=z(T)=0 



njr°°ds i 0(s i -s i _i)e- m2 ( ai - ai - 1 ) J Vz(t)exp - l -j*dtz 2 (t) 



ig 



dtz ■ A 



(6) 



1 By abuse of language the characterization 'Lorentz group' accomodates the Euclidean case as well. 
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Turning our attention to the expectation value of the Wilson exponential let us introduce 
the quantity C[z] by 

pT 



expC[z] = ( Pexp 



ig I dtz ■ A 



(7) 



J / A 



To the extent that C[z] is expected to condense the (soft) dynamics of the system, one 
should develop a strategy for its computation. To this end, we adopt the cumulant expansion, 
originally proposed in Ref [20] in connection with stochastic processes and whose significance 
for field theoretical applications has been extensively employed, see, e.g., Refs [12-16, 21]. 
Denoting cummulant correlators with double brackets, one writes 



expC[z] = exp 



w (in) 1 r 1 r 1 

E^r / dt r-" *ivW • • • z^mham ■ ■ ■ ApAti))) 



(8) 



which reduces to an identity for the case when each term in the sum finite. 

The cumulant expansion obeys a set of recursive relations in terms of field theoretical 
correlators, having the form (Lorentz and color indices for the vector field omitted) 



({A(t))) = (A(t)), 

((A( tl )A(t2))) = (P(A( tl )A(t2))) - (A(t 2 ))(A( tl )) : etc. 



(9) 



For SVM applications, it becomes important to note that factorization rules axiomatically 
operate for higher (> 3) point gluon field strength correlators. As a final note, let us remark 
that, for the perturbative case, Eq (6) corresponds to a formal relation which needs to be 
regularized by one method or other before it aquires a concrete meaning. Once this is done 
a regularization mass appears, which induces the application of appropriate renormalization 
group equations. 

The expression for the four-point Green's function by can now be summarized as follows 



G 



1 

n 

Li=4 



dsi 9{si 



-m 2 (si 



i) 



J Vz{t)S IH ... lll [z\Trce- 1 ^ 



z(0)=z(T)=0 



where 



i[z\=\ r dtz 2 {t)+iY,p l 

4 Jo i=l 



■z{ Sl )-C[z\. 



(10) 



(11) 



The quantitative analysis in this paper will be conducted with reference to the above 
generic expression, which can be seen as playing the role of an 'action functional' for a 
'particle living on the loop'. The complexity of the problem is, clearly, determined by 
the 'interaction' term C[z]. Given our intention to extend our considerations to both the 
perturbative and the non-perturbative domain of QCD the quantity C[z], which isolates all 
the dynamics in the 'particle action functional', will be organized as follows 



C[z] — C pcr t[z] + Cnon-pert^] + Cintcrff^]- 



(12) 



Our objective is to apply the particle-based casting of QCD described above by focusing 
on long distance effects. To this end we shall employ semiclassical, eikonally-based tech- 
niques. Consider, in this connection, the variation of the 'action functional' I[z] entering Eq. 
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(7) with respect to z^{t). One writes 

H)=-^W+'gM(*-.)-^ M . (13) 
Its stationary points furnish 'classical' equations of motion the integral form of which reads 

= [ dt>K,[t,t>-z«{t')\. (14) 

with the integral kernel reading as follows 

^[t,t';<(t')] = -2A(t,t')^|. (15) 

In the above equation z°{t) enters as the solution of the linear part of the 'system' and 
is given by 

*S(f)=2i£jvA(*»*)» (16) 

■i=i 

where 

A(t, t) ee -t) + ^=^V - t). (17) 

and corresponds to a straight line trajectory. 

Our expectation is that we can get a good estimate of the final result via the appplication 
of an iterative procedure on Eq (14), given the presence of a large scale, furnished by the 
incoming energy, on which an asymptotic calculation can be based. From a geometrical 
standpoint the classical solution given by Eq (14) has a profile which can be characterized 
as being of an 'eikonal type'. Our whole effort in this paper basically amounts to assessing 
first order contributions from the iterations coming both from the perturbative and the 
background (non-perturbative) gauge field sector. 

For a straight (world)line segment, joining two consecutive vertices specified, e.g, by the 
interval [sk-i, Sfc) in fig- 1 one writes (s k denotes its midpoint) 

3 3 

J2?i^H t k, si) = J^Pi^Msk, Si) = q kil (18) 
i=i i=i 

with t k G [sfc-i, Sk). Accordingly, the corresponding equations of motion read 

4 = 2iq kll (t k -8 k . 1 ) + ^(8 k . 1 ) 1 *J(0) = *J(T) = 0. (19) 

The above, 'zero' order solutions of the 'equations of motion', describe quark 'propagation' on 
straight-line (eikonal) contours. In the following section we shall use them as input solutions 
for initiating first order perturbative contributions to the four-point Green's function in the 
Sudakov kinematical region. 

3. Perturbative contributions to self-energy and vertex corrections -Sudakov 
behavior 

In this section we shall, with reference to the four-point Green's function, consider pertur- 
bative correction stemming from: (a) virtual gluon emission and absorbtion off an (eikonal) 
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fermionic line and (b) virtual gluon exchanges between such lines accross a given vertex. The 
no-recoil approximation associated with our restriction to straight line segments is equivalent 
to the 'absence' of gluon radiation. The kinematical region corresponding to large momen- 
tum transfer across a given vertex produces the conditions which give rise to the Sudakov 
form factor [1]. Going a step further, we shall determine the kinematical region for a high 
energy, four-point process which secures Sudakov conditions for each vertex. This will lead 
to Sudakov suppression behavior that is consistent with well known results [2-4]. 

3a. Perturbative self-energy and vertex corrections 

In the perturbative sector the expression for C[z], to second order of the coupling con- 
stant, is given by 

r< 2 >r*i- 1 <?c^ d t( d i\ fat fat (20) 

C pcTt [z] ---g C F ^j- 2 T - 1 j J q d tl J q dt2 lz{t2) _ z{ti)lD _ 2 . (20) 
Accordingly, the following second order expression is determined for the integral kernel: 

^i 2 Lt = 9 2 C F -^- 1 ^r2(D -2)A(t,t') fdt" . 

x[z c \t") ■ z c \t'){z*{t") - 4{t')) - z d (t") ■ At>)(z«(t") - z«(t>))}. (21) 

Generally speaking, K^ peit will induce, through an iterative procedure applied to Eq (14), 
the deviation of the classical solution from the straight line configuration. For the purposes 
of the present section, in which we shall remain within the perturbative domain of QCD, only 
leading order corrections will be considered. This, as has already been mentioned, amounts 
to employing the eikonal approximation. 

Consider the quantity 

3 



2* P 4ir D / 2 V 2 



-g 2 C F x g y n/9 r - - 1 / dt 2 dh 



Sk-l •> Sk-l 



J2 Pl (A(t 2 , Si )-A(t 1)Si ))p-2' 
i=i 



(22) 

which describes, to first perturbative order, a self energy insertion to the quark worldline 
path segment between vertices k — 1 and k. Referring to Eq (18) one finds, after some simple 
algebra, that the denominator inside the integral equals — ti). This leads to the result 

1 2 _ {4ji 2 ) 2 ~ D / 2 / D \ f s k s k -i /•«*-«*-! 



1 o [4:11 ) ' ( ' L) \ C S k~Sk-\ fSk-Sk-1 



t 2 - h 



I D-2 



= ~^C F - - ^C F \n (4ne-^VL 2 k ) + 0(e), (23) 

where 7^ is the Euler-Mascheroni constant, e = 2 — D/2 and L\ = ql(sj, — Sk-i) 2 - 

Turning our attention to gluon exchanges between adjacent quark lines, accross a given 
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vertex, we find ourselves having to deal with the quantity 

i=i 

(24) 

For subsequent purposes we also introduce the quantities 

_ q k ■ q k +\ , r\ 

Ck,k+i = I — fj 1 ■ 

|9fc||9fc+i| 

As shown in the Appendix, the I k _* k+1 can be brought into the following form: 

k+1 = - £ ^C F V 4 (c k , k+1 ) + ^C F V 4 (c k , k+1 )\n (4ne-^ +2 ^Ll) + ^C F f(c k , k+1 ) + 0(e), 

(26) 

where 2 

1 c _ c 2 i 

VAc) = — arctg = -dcotd, cos# = c (27) 

2 VI — c 2 c 2 

with 6 the angle, in Euclidean space, formed by the quark worldline paths meeting at the 
vertex k and where 

/( c ) = lim [T/ 4 _ 2£ (c)-y 4 (c)]. (28) 

Going over to Minkowski space-time entails the transcription 6 — > —27, which gives 
OcotO — > 7coth7. Denoting the total perturbative contribution from self energy (se) and 
vertex (v) corrections by P k ,Y {= h-\,k + 4-i fc fc+i + h+i fc-i)> one obtains 



rse,v 

4 "27 



CF(7fc,fc+iCoth7 fejfc+1 -l) Q + ln(47re 7B )^ + ^c F (7 fcifc+ iCoth7 fciA . +1 -l)ln^| , (29) 



where 

(30) 



-2 

A" ' 



2 - L k exp 



Consider, now, the case for which 



2 + 2- ^M+i) 



7 feifc+ iCoth7 feife+ i 



\Qk 


Qk+l\ 


\Qk\ 


\Qk+l\ 



= cosh7 fc>fc+ i» 1, (31) 

i.e. where a large momentum is imparted at the vertex. Such a situation gives rise, under 
the circumstance that gluon radiation is unobserved, to the Sudakov form factor. As already 
mentioned, the aformentioned unobservability is accounted for by the smooth (eikonal type) 
spinorial lines entering and leaving the vertex. This is consistent with the obervation that, 
in the perturbative context, the emission of a gauge field quantum from off the worldline 
contour of a matter particle induces a point of derivative discontinuity [18], whereas we have 

2 For notational economy we let c stand for the generic c k ,h+i- 
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been restricting our considerations to smooth worldline contours. In passing we note that 
the kinematical limit specified by Eq (31) leads to 



/(7fc,fc+i) - \^lk+i + \lk,k+i ln {~J~) ~ 7k,k+i> ( 32 ) 

which in turn gives 



1/2 

^2 ~ L k+1 L k exp ^7fc >fc+ i ) ^ L k+1 L k ( 2 ™* , , ) . (33) 



i n 
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9fc+i| 




9fc+i 



From the above follows that that the finite part of the vertex function, cf. Eq (29), 
assumes the form 



9fc 


Qk+i\ 


kk\ 


9fc+l 



1/2 

i 2 J^_^LLil I L k+ iL k 



\Qk 


Qk+l\ 


kk\ 





(34) 



Of course, what we are presently interested in is not the kinematical region pertaining 
to a single vertex function but the one associated with the four-point process as a whole. 
This means that the sought for kinematical region for the four-point process is one where the 
Sudakov conditions prevail for each single vertex. This is the task that we shall undertake 
in the following subsection. 

3b. Sudakov suppression for a four-point process in pQCD 

Consider the four-point process depicted in fig 1. We introduce parameters x iy i = 
1, 2, 3, 4 as follows 

Sl s 2 -si s 3 -s 2 s 4 - S3 /ocr x 

x l = —,x 2 = — , x 3 = — , x 4 = — (35) 

on the basis of which one determines (see fig 1 for designations of the momenta) 

qi = (1 - Xl )p A - x 2 (p B ' - Pb) + x 3 p B , (36) 
while the rest of the qi, % = 2, 3, 4 are specified by the relations 

92 = qi ~ PA', 93 = 9i - (Pa + Pb), 94 = 9i - Pa- (37) 

What is a priori given is that the incoming energy s = (p A + p B ) 2 sets the large scale 
for our problem. Let us then take p\ — p\ — p\, = p\, = M 2 s 2 , which means that 
Pa - Pb = Pa 1 • Pb 1 — s /2, Pa • Pa 1 = Pb • Pb' — —t/2 = |s(l — cosO) (6 the scattering angle 
in the cm frame) and pa> • Pb = Pa - Pb 1 — —u/2. With these specifications one can easily 
deduce that the conditions of Eq (31) for k — 1, 2, 3, 4 cannot be met unless three of the Xi 
are very small. Introducing a mass scale m 2 <ti s the desirable kinematical region can be 
achieved by setting x 2 — x 3 ~ x 4 ~ m 2 /s, X\ ~ 1. Simple algebraic manipulations show 
that for such a choice 

q\ ~ — , 92 ~ 93 ~ s , 94 ~ ™ 2 ( 38 ) 
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and 



gi • q 2 \ ~ m , \q 2 ■ q 3 \ ~ s, |g 3 • g 4 | ~ s, \q 4 ■ q x \ ~ m 



(39) 



from which one trivially verifies that a large momentum transfer takes place at each on of 
the four vertices. Explicitly, one deduces that 



cosh 7 12 ~ cosh 723 ~ cosh 734 ~ cosh 7 41 ~ ( — - ) 

\m 2 J 



1/2 



(40) 



With the kinematics in place, one is in position to commence with the computation of 
the amplitude. Consider, first, the contribution from the free part of the 'particle action 
functional' pertaining to the four-point contour. One obtains 

ifi=\i T dt{z»(t)y+ l j: Pl .z\s i )=i: e PiPjG^sj) 

4 J ° i=l i=lj=i+l 

S t U 

^ ^[G(s 1 , s 2 ) + G(s 3 , s 4 )] + - [G(s 2 , s 3 ) + G( Sl , s 4 )] + ~[G(s 2 , s 3 ) + G(s 2 , s 4 )], (41) 

where G(s«, Sj) = — Sj\(T — |sj — Sj\) and where we have adopted the convention pi = 

PA', P2 = Pb>, P3 - 
finds 



'Pb, Pa = ~Pa- For the kinematical region under consideration one 



4 0) ~Tm 2 . 

± free 



(42) 

This behavior can be interpreted as signifying the presence of an infrared cutoff for contours 
with length larger than 1/m and can thereby be used as the defining relation for the scale 



m 



Taking into consideration the above specifications Eq (29) assumes the following form 



2tt 



s \ 2 



rrr 



In 



A* w\\qk+i\ ( «_y 
2 



m 2 m 2 



+ 2 ln(Tm 2 ) + ln(x k x k+1 ) 



(43) 



Summing contributions from all four vertices one finally obtains 

,2 



jsc,v ^s_ c ln 
27T 



21n 



s \ 4 
"T2 



+ 41n(Tm ) + ln(xiX 2 x 3 X4) 



(44) 



Going, now, to the amplitude we note the following. First, integrations over the param- 
eters x 2 , x 3 , X4 in the region x 2 ~ x% ~ a: 4 ~ m 2 /s should be performed and second, the 
computed (up to first perturbative order) contributions to the 'particle action functional' 
enter as exponentials. We proceed with their assessment by considering one by one the three 
terms composing J se ' v . 

The contribution coming from the last term of the above equation gives rise to integrals 
of the form 



y o dx k x k 



1 (J. 

l-f^ln(^) Vm 



-{i-S^m(^r)} 



(45) 
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ensuring that the kinematical region under consideration does not have zero measure. More 
specifically, as will be argued in the next section, the strong coupling constant can be defined 
at a scale for which 

^C F 0(1), (46) 

hence it should eventually give non-negligible contributions. 
The second term gives 



cxp 



^C F ln (-L 41n(Tm 2 ) 
2-7T \m 2 J 



(Tm 2 ) 



2\-4^C F \n- 

Zty t 



(47) 



Integration over T, as entailed by the worldline expression, leads to a final contribution of 
the form 



poo aQ 

/ dTT 3 (Tm 2 y 4 ^ 
Jo 



2x-4^C F ln^ 3r -Tm 



/ °rfT(T)- 1+4 ( 1 --^ ln ^)e- Tm2 
Jo 



4(1- ^C F \n -4 



2tt 



(48) 



which is of finite order. 

The notable contribution at the perturbative level is furnished by the first term, i.e. 



G ~ cxp \ -—C F \n (A, ) In 



7T 



i n 

S \ 4 



m 2 \m 2 



(49) 



and displays Sudakov suppression for an exclusive, high energy four point process. 

As a final note, let us observe that the presence of the regularization mass scale \x implies 
an eventual renormalization group running, given the final result should not depend on it. 
Such a running calls for 'initial conditions' set at some minimum scale \i m %n- We shall defer 
the discussion of this matter until non-perturbative deformations to the integral kernel will 
have also been taken into account. 

4. Non-perturbative considerations: Deformation of the straight line contours 

As mentioned in the introduction, non-perturbative effects associated with the four-point 
Green's function under study will be assessed through the employment of the Stochastic 
Vacuum Model (SVM) [12,13]. Its central objective is to incorporate established observa- 
tions/results regarding the structure of the QCD vacuum whose starting point can be traced 
to Ref [21]. The SVM attempts to summarize all that is known and/or surmized about the 
properties of the QCD vacuum through a set of three axioms which are expressed in terms 
of field strength, as opposed to field potential, correlators. The underlying stochasticity 
assumption for the vacuum state facilitates the direct application of the cumulant expansion 
-see, e.g., relevant review articles [22,16,20]. A concrete, as well as practical, way to apply 
the SVM scheme to specific situations is by using the background gauge fixing method [23], 
with the background gauge fields becoming the agents of the non-perturbative dynamics. 
Specifically, one employs the gauge potential splitting A® = a® + where the a® are 
associated with the usual perturbative field modes. The B®, on the other hand, enter as 
dynamical fields, assigned with the task of carrying the non-perturbative physics through 
field strength correlators which obey factorization rules according to which higher order 
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gluon field strength correlators are expressible in terms of two-point ones (third axiom of the 
SVM). 

Some preliminary matters should be dealt with from the outset. First, let us mention 
that we shall be employing throughout the Fock-Schwinger (F-S) gauge [24,7], namely 



rx rl 

B a( x ) = ~ du v {d^u p )Fp U {u) = -(x - x ) p j da aF^ u (x + a(x - x )), (50) 

Jx JO 



which facilitates the passage from gauge field potential to field strength correlators. Specifi- 
cally, one has 

(gB tl2 (x(t 2 ))gB lll (x(t 1 )))B = (x2-xo) V2 (x 1 -xo) l , 1 [ da 2 a 2 [ da x a x 

Jo Jo 

x (gF^ 2U2 (x + a 2 (x 2 - x ))gF^ lul (x + a^X! -x ))) B . (51) 

Upon setting Ui = x + aix(t), i — 1, 2 one determines 

^rc(gF; 2V2 {u 2 )gF; ivi { Ul )) B = 2N C 'A^ m (u 2 - u x ) 

= 2tr c ((f)(x , u 2 )gF^ 2U2 (u 2 )<f)(u 2 , x )(f)(x , u^gF^ (ui)0(ui, x )) B , (52) 



where <J)(xo,Ui) = Pexp [ig f*° dv ■ B(v)j. In the F-S gauge this factor is unity. Its insertion 
serves to underline the gauge invariance of the field strength correlator. 

Consider, now, the second term entering the rhs of Eq (12). Substituting the label 'non- 
pert' by 'bkgd', given that the non-perturbative dynamics of QCD are carried by the gauge 
field modes B° we write for the main object of computational interest 

Cbkgd = ~ Cf N 2 Z i J da ^ a 2 J o dt 2 z(t 2 ) da! a! dhz^A^^ [a 2 z(t 2 ) - anzfa)] 



having shifted the reference point x to zero, so that Ui = aiiz(ti), i — 1, 2. 

A convenient representation of the correlator, to which we shall be referring throughout 
our analysis, is [22,16,20] 



+ 



1 d 



2 dzp il ... 



1 d 



ZlA2^V-2V\ Zv-2$ U,-2V\) D \[Z ) + (^2^M2Ml %112$W2IJ,\)D\(^Z ) . (54) 

Z oz v 



Now, the integral kernel, in terms of which the lowest order non-perturbative corrections to 
the amplitude will be determined, is given by 

* sd (0] = "2A(t, t')jj^f$jz d } (55) 

so, naturally, we need to focus our attention to C^ kgd [2;]. 

Given the large energy scale s entering our problem, we can safely assume that any ratio 
formed between some other (energy) scale characterising the system intrinsically, such as the 
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string tension or the inverse correlation length (squared) and s will be very small. It can be 
demonstrated that, under these circustances 



^1 



bkgd 



[A 



<j S m - U 



=-- r 

2 Jo 



z ■ z) 2 — Z 2 Z 2 , 



(56) 



where a is the string tension, entering, according to the premises of the SVM, through the 
relation a = \ J d 2 zD(z 2 ), while S min rpresents the minimal surface bounded by the Wilson 
loop associated with the closed quark contour. For the particular configuration given by 
the straight line segment solution z®, with reference to our parametrization of the four-point 
function loop, one obtains 



CbkgdN « -2a{si(s 2 - Si)\qi ■ q 2 \ + (s 2 - s 1 )(s 3 - s 2 )\q 2 • ^3 1 } 



(57) 



from which one surmizes that the corresponding corrections to the linear (eikonal) result will 
be ~ a I s. 

Referring, now. to Eq (14) the first order non-perturbative contribution to the 'particle 
action' reads 



4{t) = zl{t) + 2a jf dt'A(t,t') S -^^ + O{o\g\ 



One determines 



{+>\/ -O/.tx^i 



K(t>) 



where 



<7<°) = (z° • i ) 2 - (z°) 2 (z ) 2 . 



(58) 



(59) 



(60) 



It is a straightfoward matter to deduce that, for the kinematical region under consider- 
ation, the the intervals (0, s±\ and (S3, S4] remain linear while the other two segments are 
deformed. Specifically, one obtains 



and 



t'e(0, Sl }: z c J(t) = 2iq llx t, 

t' E (si, s 2 ] : z*(t) = 2iq 2il {t - Sl ) - 2iq 2fl a(t - Si) 2 + z£(si), 
t' E (s 2 , s 3 ] : z*(t) = 2iq 3>1 {t - s 2 ) - 2iq 3fl a(t - s 2 f + z*(s 2 ), 

f G(s 3 ,s 4 ]: z*(t) = -2iq 4fl (T-t). 



(61) 

(62) 
(63) 

(64) 



In the above relations the following auxiliary 'momentum' variables have been introduced 



q2n = (hii 



gg 

<?2 ■ <?3 



5i • I2 



s 2 — s l 



(65) 
(66) 
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while the re-adjusted momenum variables are given by 



qip = qin + 2a 



<?2 M s 2 - si - 



IT 



Q2n = <?2 M + 2a 

9V = qip + 2a 



92 M \S2-Sx 



£2 f? 

2T 



+ 93/. S3 - S 2 ~ 



+ 93/. S3 - S 2 ~ 



S 3 s 



2T 



s 2 - s 2 " 

*3 S 2 



2T 



-?2 M " 



and 



94 M = 94/< + 2(7 



2T 



"92/, 



+ 93/. S3 - S 2 



S3 - s 



2T 



s 2 - s 2 
2T 



93/. 



s 2 - s 2 

a 3 ^2 

2T 



(67) 

(68) 
(69) 

(70) 



The resulting worldline configuration for the four-point (Sudakov) process is depicted in fig 
2. 

We now proceed to substitute the above 'deformed' line paths into the second order 
expression for Cl^Jz] given by Eq (20) and asses the new worldline configuration associated 
with the four-point process (for Sudakov kinematics). The expectation is that the non- 
linear deformations will produce corrections of the order of cr/s, while the linear ones will 
simply adjust the corresponding expressions to the displaced momentum variables. For the 
kinematical region under consideration and for a < m 2 the results of section 3 remain intact, 
while for a > m 2 we determine 



q 1 ~ a- 



m 



q\ ~ cr, ql ~ s, q\ ~ m 2 



and 



|gi-g 2 |~a, \q 2 ■ q 3 \ ~ s, \q 3 ■ q 4 \ ~ s, \q 4 ■ q±\ ~ a. 
From the above, the following estimates for the quantities cosh 7^+1 can be deduced 



(71) 
(72) 



\qi 


h\ 


\qi\ 


\h\ 


\os 


94 1 


\h\ 


|94| 



s 



\h 


93 1 


\<ht\ 


\h\ 


I94 


9i| 


l&l 


\qi\ 



a 



(73) 



Having adjusted our expressions to the new situation, which has brought into our analysis 
input from the non-perturbative (confining) structure of the QCD vacuum, we are ready to 
sum the soft effects associated with the four-point Sudakov process. To this end we shall 
turn to the renormalization group equation, which will facilitate the running from a low, 
'infrared' scale /i min up to the largest energy scale implicated for the process, namely s. This 
task will be carried out in the section that follows. 

5. Renormalization group running 

Repeating the analysis of section 3, one can use the results of the previous subsection to 
obtain the modified expressions for the J^ e ' v , k — 1, ...,4. We present the final result for the 
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four-point process, replacing the one given by Eq (44), which reads as follows 



jse,v _ ^lc F ln 
2tt 



2 In 



m^fa \m 2 



+ Aln(Tm\fa) + \n(x 1X2X3X4) 



2n 4 



a 



(74) 



The last two terms in the curly brackets are connected with T- and x^-integrations, just as 
in the perturbative case. We now introduce the mass scale Umin 



d 



(9 In s 



m 



o \ 1 

s. 



(75) 



to designate the point at which initial conditions for the renormalization group running are 
to be set. As noted in the beginning of the present section it should, of course, hold that 
/"miii > A-qcd- The limit a — > m 2 gives the perturbative estimate /z min = jfe ■ 

With these specifications let us proceed with renormalization group considerations. Re- 
ferring to Eq. (6), which gives the unrenormalized expression for the four-point function, we 
write, in condensed notation, 

[ M, KU] = J Vz(t)S lit ... ltl [z]exp ~ £dtz 2 (t) -iY,Pi- 



z(Q)=z(T)=Q 



x ( Tr c Pexp 



ig dtx-A) = Z[{pi},{si},/i,e]^2. m [fe},{si},^(/i), //], (76) 



where labels (0) and (R) denote unrenormalized and renormalized quantities, respectively. 
We find it convenient to run the quantity 



d 



dins 



dins 



dins 



(77) 



from /i min to s. The renormalization group equation, which reflects \i independence, reads 



1 d i nU (R) = _ dd \ nZ 



-2T 



dln/j Ins dln/i Ins 

The anomalous dimension can be read off Eq (82): 

T = ^C F + 0(a 2 s ). 

Retaining leading contributions we present the solution in the form 



(78) 



(79) 



t 



hu7(«V = s) = -f jln (£) r[a s (t)}-2 J™ f In ( -f- ] T[a s (t)] + lnU^(^ = 

ft min \r^miil/ 



(80) 
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which leads to 



4C 



s = 



F 



A) 



in 



In I 



\n 7 - 



In 



A 2 



i + 



In 



U 2 . 
mm 

A 2 



e-in 



A 2 



+lnf/^(/i 2 



From Eq (74) one determines 



2 \ 
A* mini 



(81) 



2 

A^min 



2 - ) 

mm/ 



In ( — ) [4 ln(Tmv/ff) + lnfcia^au)] + sbldg, (82) 

7TZ / 



where 'sbldg' stands for subleading terms. 

If one were to go back to our full expressions and perform the integrations over the 
Xk and T, one would realize that the coupling as — as(/-f mm ) * s consistent with the claim 
-, as stated by Eq (46). Accordingly, the nonperturbative contribution re- 



■Z^F i n ( 

adjusts the Sudakov result for the four-point process as follows 



A'min 



4tt 1 
I3 ln(s/A 2 ) 

relation. It would then follow that 




+ sbldg. 



(83) 



Suppose one sets a,s{s) 



and keeps terms to first order in as in the above 



Xj{R) 



cxp 



Os_ 

'2?r' 



In 2 



In 2 



(7 



(84) 



a result which explicitly displays a visible modification brought about through the inclusion 
of non-perturbative input to the analysis of the four-point process. On the other hand, 
it is observed that this modification is not strong enough to negate the role that Sudakov 
behavior plays in the analysis of exclusive processes in QCD. 

6. Concluding remarks 

In this paper we have defined the Sudakov, QCD mediated, four-point process as one for 
which: a) A large momentum transfer takes place at each single vertex and b) no observable 
gluon radiation, in the form of jets, occurs above a given infrared scale. Such an 'idealized' 
setting allows for a direct study of the Sudakov suppression in QCD in a way which is 
liberated from extra phenomenological burden. The key to atttaining our goal was to isolate 
that corner of the kinematical space which realizes the Sudakov conditions for the process 
as a whole. Granted, in a realistic situation, such as one pertaining to meson-meson elastic 
scattering, the corresponding kinematical analysis is far more complicated. On the other 
hand, our 'idealized' setting is in position to capture the essence of the Sudakov kinematical 
constraints and lead to results which check directly with the ones corresponding to the more 
'realistic' analysis. In short, our claim is that the corner of the kinematical space identified 
in this work corresponds to the one which defines the Feynman picture, as opposed to the 
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quark counting rule, for hadron-hadron elastic scattering. The clear dividents of the present 
approach is that it enabled us to get a glimpse on non-perturbative alterations to the Sudakov 
result, as one approaches Aqcd from above, through the employment of the SVM. 

The following comment, given in retrospect, might be useful. It concerns the utilization of 
the scale m 2 for isolating the relevant integration region for Sudakov kinematics, cf. Eq (42). 
The underlying hypothesis is that in a complete computation, i.e. with the nonperturbative 
contributions fully taken into account, such a scale would be determinable in terms of the 
dimensional parameters entering the problem: m 2 = m 2 (M 2 ,a, ■ • •). The nature of the 
present calculation was such that it renders the non-perturbative corrections in the region 
m 2 > a negligible, whereas they become more important for m 2 < a, cf Eq (84). The point 
is that our present effort has contained itself to the perturbative domain of the theory, albeit 
bringing in some non-perturbative input (interference) through the SVM. We anticipate that 
the extension of the present analysis to the Regge kinamatical region (always for the four- 
point function) will present richer features due to the fact that the relevant gluon exchanges 
will occur among opposite, so to say, (world)line segments, i.e. segments separated by two 
consecutive vertices. Our next goal is, indeed, to extend the analysis of this work, for the 
same process, to the Regge kinematical regime. 
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Appendix 

In this Appendix we derive Eqs (26)-(28) in the text, starting from Eq. (24). The latter 
becomes, via the use of Eq (18), 



T K 
1 k-i 



k 



, +J - - 2 9 c F —^^r - ij qk+rqk J dtl J dt 2 ~ 



ID— 2 



Eft (A(t 2 , Si ) - A(ti, Si )) 

i=l 

(A.l) 

The denominator of the expression inside the integral can be simplified, once appealing to 
Eq (18): 

3 

J^Pv(A(*2,s*) _ A(ti, s^) = qk+iA^ +Pfc M ( s fe ( A - 2 ) 
i=i 

for t 2 e [sfe+i, s fe ), ii G [s fc , s fe _i) 
One thereby obtains 

(A.3) 

Reparametrizing according to s k — ti — > t[, t 2 — s k — > i 2 an d noticing that g^+i +Vk = % 
one finds 



D-2 



- ? 2 ^ r (f - ^^"^^rb (A - 4) 

with the Cfc^+i as defined in Eq (25), to be written simply as c from hereon. 
In the above result we have introduced 

V D = ^ D (x, c) + Ix 4 - D A D (i, c), (A.5) 

where 

Lk+l \Qk+l \\ s k+i — s k\ / a m 

x = ~~ ; — = — i — n r (A.b) 



and 

An(* d = C 2 F f 1 D - 2 . D ~ 1 . 1 C A X + ° F ( X D -^ D -^. 1 ~ C2 * 

(A.7) 

For D = 4 one obtains 

c ^/ 1 — c 2 c ^/ 1 — c 2 
c) = , arctan . arctan , (A. 8) 

which, upon substitution in Eq (A5), leads to the desired result. 
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Figure captions 



Fig. 1. The four-point process under cosideration in the eikonal approximation (straight 
worldline segments). 

Fig. 2. Deformed four-point worldline contour on the account of SVM induced non- 
perturbative effects (still, slightly above Aq C d)- 



20 




1 




1 



